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Abstract. From the analysis of the C + fine-structure population ratio in the damped Lya system at z a b s = 3.025 
toward the quasar Q0347-3819 we derive an upper bound of 14.6 ± 0.2 K on the cosmic microwave background 
temperature (Tcmbr ) regardless the presence of other different excitation mechanisms. The analysis of the ground 
state rotational level populations of H2 detected in the system reveals a Galactic-type UV radiation field ruling 
out UV pumping as an important excitation mechanism for C + . The low dust content estimated from the Cr/Zn 
ratio indicates that the IR dust emission can also be neglected. When the collisional excitation is considered, we 
measure a temperature for the cosmic background radiation of Tcmbr =12.lt_\'2 K. The results are in agreement 
with the Tcmbr = 10.968 ±0.004 K predicted by the hot Big Bang cosmology at z a ^ s = 3.025. 
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1. Introduction 

In the standard Big Bang model (SBB) the tempera- 
ture of the relic radiation from the hot phase of the 
Universe is predicted to increase linearly with redshift: 
Tcmbr {z) = T CM br (0) (1+z) (e.g., Peebles 1993). At the 
present epoch direct measurements show that Tcmbr (0) 
= 2.725±0.001 K (la c.L), and that the relic radiation fol- 
lows a Planck spectrum with very high precision (Mather 
et al. 1999). 

As pointed out by Bachall & Wolf (1968) the CMBR 
temperatures at earlier epochs can be measured from the 
analysis of quasar absorption line spectra which show 
atomic and/or ionic fine-structure levels excited by the 
photo-absorption of the CMBR. Among the species with 
fine structure levels the C 1 and C II show an energy sepa- 
ration, from 23.6 K up to 91.3 K, which make them sensi- 
tive to the CMBR, in particular as the redshift increases. 
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However, C 1 is generally fully ionized and rarely detected, 
while the C 11 ground-state transitions are strongly sat- 
urated, thus making column densities rather uncertain. 
In addition, non cosmological sources (such as particle 
collisions, pumping by UV radiation, IR dust emission 
and by other sources) may compete with the CMBR to 
populate the excited fine-structure levels. Only indepen- 
dent knowledge of ambient radiation field and of parti- 
cle densities allows to disentangle the contribution of the 
background radiation from that of other mechanisms. For 
these reasons previous measurements place upper limits to 
Tcmbr rather than real measurements, albeit quite strin- 
gent ones (Meyer et al. 1986; Songaila et al. 1994; Lu et 
al. 1996; Ge, Bechtold & Black 1997; Roth & Bauer 1999; 
Ge, Betchold & Kulkarni 2001). 

Recently, Srianand, Petitjean & Ledoux (2000) from 
the H2 analysis in the DLA at z a b s = 2.3371 toward the 
quasar Q1232+0815 were able to infer the UV radiation 
field in the absorber. Then by means of C I, C I*, C I** , C II 
and C 11* they obtained a Tcmbr = 10 ± 4 K, while SBB 
predicts Tcmbr = 9.09 K. However, the H2 abundance 
measurement at z a b s = 2.3371 by Srianand et al. (2000) is 
in contradiction with their recent estimation of the deuter- 
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Fig. 1. The Galactic interstellar radiation fields at galac- 
tocentric distances Dq — 5 and 10 kpc from Mathis et al. 
(1983) are compared with the black body spectra calcu- 
lated at different redshifts using the relation Tcmbr (z) 
= Tcmbr (0) (1 + z). The positions of the excited levels 
of C° and C + , suitable to constrain the Tcmbr values at 
different redshifts, are indicated by vertical dashed lines 

ated molecular hydrogen abundance (Varshalovich et al. 
2001). The ratio HD/H 2 ~ (0.8 3.0) x 10~ 3 , whereas it 
is 10~ 7 10~ 6 in the ISM diffuse clouds (e.g. Wright & 
Morton 1979). Until this discrepancy is clarified the value 
of Srianand et al. (2000) should be taken as an upper 
limit of Tcmbr < 14 K at 

z abs — 2.3371. So far all mea- 
surements have been found to be consistent with the SBB 
model prediction. 

In this letter, we present a new measurement 
of Tcmbr at higher redshift, z a b s = 3.025, from the 
VLT/UVES spectra of Q0347-3819. 

2. Data analysis and results 

The spectroscopic observations of Q0347-3819 obtained 
during UVES commissioning at the VLT 8.2 m telescope 
are described in detail by D'Odorico, Dessauges-Zavadsky 
& Molaro (2001) and by Levshakov et al. 2002 (LDDM, 
hereinafter) . 

In LDDM relevant physical properties for the damped 
Lya system (DLA) at z a b s = 3.025 are obtained by ana- 
lyzing numerous H 2 and metal absorption lines associated 
to the DLA. The z a b s = 3.025 system exhibits a multi- 
component velocity structure spanning over ~ 80 km s _1 . 
The main component at z a b s = 3.024855 has a hydrogen 
column density of TV (Hi) = 2.52(±0.04) x 10 20 cm" 2 and 
shows the presence of molecular hydrogen with a fractional 
abundance of /(H 2 ) = 3.25(±0.17) x 10~ 6 . Several neu- 
tral and ionized species associated with this cloud have 



been analyzed. In particular the Cn 1036.3367 and Cll* 
1037.0182 lines have been identified. 

A line absorption model for the z a b s = 3.025 system 
which is able to reproduce the line profiles for the whole 
set of atomic and ionic species has been elaborated in 
LDDM. This model allows to get a reliable column den- 
sity of the saturated lines. LDDM obtained a column den- 
sity of iV(Cn) = 5.05(±0.28) x 10 15 cm" 2 for the Cll 
1036.3 main component. This C column density is con- 
sistent with what can be inferred from the other elements 
measured in the system by means of unsaturated lines. For 
instance, if C goes in lockstep with the undepleted Znii 
we would obtain iV(Cn)= 3.1 xlO 15 cm~ 2 , assuming so- 
lar photosphcric values from Grevesse and Sauval (1998), 
while we would obtain 7V(Cn)= 5.8 x 10 15 cm~ 2 if C fol- 
lows Ar I, with the Ar solar value quoted in Sofia & Jenkins 

(1998) . Neutral carbon is not detected and N(C i)/N(C n) 
< 7.9 x 10~ 5 . 

The column density for the iV(Cll*) 1037.0182 main 
component is 2.26(±0.12) x 10 13 cur 2 . Prochaska & Wolfe 

(1999) reported the detection of the Cll 1334.5323 and 
Cll* 1335.7077 lines in the same system. For the latter 
line, which is unsaturated, they provide a column density 
of AT(Cll*) = 3.00(±0.23) x 10 13 cur 2 , which refers to the 
total system. When we correct this value according to the 
relative ratios between the components [1 : 0.195 : 0.044 : 
1.952 (LDDM)], we obtain for the main one iV(Cll*) = 
1.83(±0.14) x 10 13 cm" 2 . The C II* 1335.7077 line is likely 
blended with the Cn* 1335.6627 which produces the blue 
asymmetry present in the Keck spectrum at -10 km s -1 
(cf. Fig. 5 in Prochaska & Wolfe). The relative strengths 
of the two blended transitions is f 1335.1 1 'J '1335.6 = 8.7. If 
we correct iV(Cii*) by the corresponding factor, for opti- 
cally thin lines we obtain JV(Cll*) = 1.64(±0.11) x 10 13 
cm -2 (main component). The weighted mean between the 
VLT and Keck quantities is N(C 11*) = 1.92(±0.08) x 10 13 
cm -2 . Combining this value with the ground level col- 
umn density obtained from the VLT we derive a ratio 
N(Cu*)/N(Cu) = 3.8(±0.3) x 10~ 3 . 

The ground state of the C + ion consists of two levels 
2s 2 2p 2 P"/ 2 3/ 2 with an energy separation of AE = 63.42 
cm -1 which corresponds to A — 157.7 fim. The excited 
level can be populated by several mechanisms such as colli- 
sions, fluorescence or IR photon absorption, which include 
also the CMBR. In the following we use an effective tem- 
perature T e g to characterize at A = 157.7 [im the proper 
spectral energy density of the local IR field approximated 
by a Planck spectrum with T = T e g. 

In equilibrium, the population ratio of the upper level 
n 2 to the lower level n\ , in ions with a doublet fine struc- 
ture in the ground state, is given by: 

rjv _ Qi,2 + wig ,y 
n-i Q2.1 + A 2 ,i + w 2 ,i 

where Wi^, w 2 ,i and Qi,2, Q2,i are the photo-absorption, 
radiative decay, collisional excitation and de-excitation 
rates, respectively. A 2 .i is the radiative transition prob- 
ability which is A 2 ,i = 2.291 x 10~ 6 s" 1 for C+. 
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If only the background radiation contributes to the 
population of the excited fine-structure states, eq. (1) 
gives: 



Tcmbr {z) — 



91.325 



In [g* N{CII)/gN{CII*)] 



(2) 



where g* and g are the statistical weights of the corre- 
sponding levels. Thus the estimated 7V(C n*)/iV(C n) ra- 
tio would yield T c g = 14.6±0.2 K, which can be considered 
as a firm upper limit to the Tcmbr • The value expected 
from the standard model at z a b s = 3.025 is of Tcmbr = 
10.968 ±0.004 K. Thus the observed 7V(Cll*)/7V(Cll) ra- 
tio provides a stringent limit to the Tcmbr at z a b s = 3.025 
with only little room left for extra contributions. 

In the following we show how the detection of H 2 in 
the same component where C II* and C II are observed can 
provide additional information on the presence of other ex- 
citation processes. In this discussion we assume that the 
molecular and ionic species trace the same material as it 
is suggested by the similar broadening shown by the line 
profiles and by the absence of any evidence for an associ- 
ated dense H n region gas on the line of sight as argued in 
LDDM from the non detection of the N n* 1084.580 and 
1084.562 lines. 

The H 2 populations over the J = to J = 5 rotational 
levels of the ground electronic-vibrational state provide 
an excitation temperature of T cx — 825 ± 110 K and the 
kinetic temperature is also estimated to be Tkin ^ 430 
K (LDDM). The population ratios of the higher J levels 
N(5)/N(S) and N(A)/N(2) are sensitive to the UV pump- 
ing. The measured rate of photo-absorption /3o ~ 2 x 10~ 9 
s _1 is very close to the average interstellar radiation field 
in the Galaxy. With this constraint on the UV flux the 
fluorescent excitation process has a rate ~ 9.3 x 10~ n 
s -1 , which is rather low and can be neglected according 
to Silva & Viegas (2001). 

The rates of the radiative processes u>i i2 and u>2,i may 
be caused by the cosmic microwave background radiation 
at z a bs = 3.025, but also by local sources of infrared radia- 
tion like diffuse emission from dust heated by OB stars to 
temperatures T ~ 10-^20 K as observed in the Milky Way 
(Mathis, Mezger & Panagia 1983). The possible contribu- 
tion from the heated dust is illustrated in Fig. 1 where 
the MW interstellar radiation fields at galactocentric dis- 
tances Dq = 5 and 10 kpc are shown along with the black 
body spectra calculated at different redshifts using the lin- 
ear relation T CM br (z) = T CM br (0) (1 + z). This is repre- 
sentative of our system since, as we have discussed above, 
the intensity of the UV field in the z a bs =3.025 cloud is 
found to be of the same order of magnitude as in the MW. 
The positions of the excited levels of C° and C + , which 
are suitable to restrict Tcmbr , are also indicated by ver- 
tical lines. Fig. 1 shows that the diffuse FIR reemission 
of stellar radiation by dust grains, if the dust emissivity 
at A = 157.7 yum is equal to the highest value measured 
at Dq = 5 kpc in the MW, always remains lower than 
the expected CMBR. The corresponding photo-absorption 
rate is wf 2 st — 1-8 x 10~ n s _1 , but the expected rate in- 




Fig. 2. Monte Carlo simulations of the probability den- 
sity function of T cS for the value of AT(Cll*)/iV(Cll) = 
3.8 x 10~ 3 . The mean value is T off = 12.1+^ K [the 
\{l—p) and \ {1 +p) quantiles were used to estimate the 
uncertainty interval at p — 0.95]. The Tcmbr from the 
standard Big Bang cosmological model is marked with a 
vertical dashed line 



duced by the relic radiation is Wi2 ltlti — 1.1 x 10~ 9 s _1 . 
Moreover, LDDM estimated from the [Cr/Zn] abundance 
ratio that the dust content in the z a b s = 3.025 absorbing 
region is about 30 times lower as compared with the MW 
mean value, so that we may exclude significant contribu- 
tion from dust emission. 

We now consider the information on the particle den- 
sity, since the upper level of C + may be populated by 
collisions with several particles such as electrons, e~, hy- 
drogen atoms, H°, protons, H + , and molecules, H2 . The 
J=2 level of H2 has a rather long radiative lifetime and 
is the more sensitive to the collisional de-excitation. The 
critical density above which collisional de-excitation be- 
comes important is riy=3.6 cm -3 (LDDM) and, there- 
fore, jih > 4 cm~ 3 is required to maintain the observed 
N(2)/N(0) ratio at T kin ~ 400 K. Arguments based on 
the production rate of H2 imply that the volumetric gas 
density, jih, ranges between 4 and 14 cm~ 3 . 

The H°-C+ collisional rate is of - 1A5 x 10~ 9 nH 
s _1 in the range 10 2 K < T kin < 10 3 K (Launay & Roueff 
1977). Collisions with electrons have the highest rates but 
the electron density is rather low. Electrons in H 1 re- 
gions come mainly from carbon photo-ionization, so that 
n c ~ (C/H) n H , which is ~ 3 x 10~ 5 n H for the z abs = 3.025 
system. The rate is Q\ 2 — 2 x 10~ 7 n o s _1 and therefore 



the collisional rate becomes Q\ 2 — 6 x 10~ 12 nn s _1 , 
which is much lower than that of the hydrogen collisions 
for the same temperature interval (Silva & Viegas 2001). 
H 2 molecules do not contribute to collisions considering 
the low fractional abundance measured in the system. The 
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corresponding de-excitation rate, calculated from the prin- 
ciple of detailed balance, is ~ 9.11 x 10~ 10 rtH s" 1 , 
for T kin ~ 400 K. 

We calculated the probability density function of T e g 
using statistical Monte Carlo simulations which suggest 
that the errors are normally distributed around the mean 
value of iV(C u*)/N(C n) with the dispersion equal to the 
probable error of this ratio, while jih is evenly distributed 
between 4 cm -3 and 14 cm -3 . The result is presented in 
Fig. 2. The most probable value of T e g obtained in this 
analysis is T e g = 12. l+g \ K. The lower and upper errors 
of T ff correspond to the |(1 — p) and |(1 +p) quantiles, 
respectively (the central 100p% confidence interval was 
used with p — 0.95). 

Since we have considered collisions and excluded 
fluorescence and dust emission as significant processes 
in the population of the excited levels, T e g is actually 
Tcmbr f° r this particular DLA. Thus our measurement 
of N(Cu*)/N(Cu) = 3.8(±0.3) x 10~ 3 leads to the most 
probable value of T c g = 12.1 K which is only 1.1 K higher 
with respect to the predicted Tcmbr and fully consistent 
within errors. 

In Fig. 3 all the previous estimations of Tcmbr are 
shown. Our result, together with upper limits presented 
in Fig. 3 support the linear evolution of the CMBR within 
the framework of the SBB model. 

Alternative non-adiabatic cosmological models in 
which photon creation takes place as the Universe ex- 
pands predict a different temperature-redshift relation of 
the type T CM br (*) = T CM br (0) (l + z)^) (Lima, Silva 
& Viegas 2000). At high redshift the deviation becomes 
more pronounced and our measurement set a limit to j3 < 
0.22 (2 a). 

3. Conclusions 

The analysis of the H2 lines in the damped Lya absorber 
at z aD .s = 3.025 toward QSO 0347-3819 allows us to es- 
timate the local excitation mechanisms which populate 
the fine-structure levels together with the Tcmbr ■ From 
the N(C ii*)/N(Cii) ratio we measure the temperature of 
the local background radiation of Tcmbr = 12.1 jig" \ K 
which is consistent with the temperature of the cosmic 
background microwave radiation of 10.968 K predicted 
by the standard Big Bang cosmology at the redshift of 
the absorber. 
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